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. ABSTRACT 

We present evidence that high-velocity cloud (HVC) Complex C is a low-metallicity 

. _ gas cloud that is plunging toward the disk and beginning to interact with the ambient 

' gas that surrounds the Milky Way. This evidence begins with a new high-resolution 

■ (7 km s^^ FWHM) echelle spectrum of 3C 351 obtained with the Space Telescope 

^ ■ Imaging Spectrograph (STIS). 3C 351 lies behind the low-latitude edge of Complex 

. C, and the new spectrum provides accurate measurements of O I, Si II, Al II, Fe II, 

. and Si III absorption lines at the velocity of Complex C; N I, S II, Si IV, and C IV 

r~| ! are not detected at 3a significance in Complex C proper. However, Si IV and C IV 

Oh| as well as O I, Al II, Si II and Si III absorption lines are clearly present at somewhat 

O ' higher velocities associated with a "high-velocity ridge" (HVR) of 21cm emission. This 

' high-velocity ridge has a similar morphology to, and is roughly centered on. Complex 
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C proper. The similarity of the absorption lines ratios in the HVR and Complex C 
suggest that these structures are intimately related. In Complex C proper we find 
[0/H] = -0.76+ 

a2i- ^'^^ other species, the measured column densities indicate that 
ionization corrections are important. We use collisional and photoionization models 
to derive ionization corrections; in both models we find that the overall metallicity 
Z = 0.1 — 0.3^0 in Complex C proper, but nitrogen must be under abundant. The iron 
abundance indicates that the Complex C contains very little dust. The size and density 
implied by the ionization models indicate that the absorbing gas is not gravitationally 
confined. The gas could be pressure-confined by an external medium, but alternatively 
wc may be viewing the leading edge of the HVC, which is ablating and dissipating as it 
plunges into the Milky Way. O VI column densities observed with FUSE toward nine 
QSOs/AGNs behind Complex C support this conclusion: N{0 VI) is highest near 3C 
351, and the O VI/H I ratio increases substantially with decreasing latitude, suggesting 
that the lower-latitude portion of the cloud is interacting more vigorously with the 
Galaxy. The other sight lines through Complex C show some dispersion in metallicity, 
but with the current uncertainties, the measurements are consistent with a constant 
metallicity throughout the HVC. However, all of the Complex C sight lines require 
significant nitrogen under abundances. Finally, we compare the 3C 351 data to high- 
resolution STIS observations of the nearby QSO H1821-I-643 to search for evidence of 
outflowing Galactic fountain gas that could be mixing with Complex C. Wc find that 
the inter mediate- velocity gas detected toward 3C 351 and H1821-I-643 has a higher 
metallicity and may well be a fountain/chimney outflow from the Perseus spiral arm. 
However, the results for the higher-velocity gas are inconclusive: the HVC detected 
toward H1821-I-643 near the velocity of Complex C could have a similar metallicity 
to the 3C 351 gas, or it could have a significantly higher Z, depending on the poorly 
constrained ionization correction. 

Subject headings: Galaxy: abundances — Galaxy: halo — ISM: abundances — 
ISM: clouds — quasars: individual (3C 351, H1821+643) 

1. Introduction 

The Galactic high- velocity clouds (HVCs), gas clouds detected via 21 cm emission or 
UV/optical absorption at velocities that deviate substantially from normal Galactic rota- 
tion (Wakker & van Woerden 1997), have potentially important implications regarding the 
structure and evolution of the Galaxy and Local Group. For example, if the HVCs provide 
a sufficient quantity of infalling low-metallicity gas, they can alleviate the G-dwarf problem, 
the long-standing discrepancy between the observed metallicity distribution of G-dwarf stars 
and theoretical expectations (e.g., Larson 1972; Tosi 1988; Wakker et al. 1999; Gibson et 
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al. 2002). The HVCs may have important cosmological imphcations as well. Hierarchical 
models of galaxy formation within the cold dark matter framework predict many more dwarf 
satellite galaxies within the Local Group than have been detected/identified (e.g., Klypin et 
al. 1999; Moore ct al. 1999). A variety of solutions to this problem have been proposed, 
including the possibility that the HVCs are the missing dwarfs which, for some reason, have 
not formed readily detectable stars (e.g.. Blitz et al. 1999; Klypin et al. 1999; Gibson et al. 

2002) . 

However, the nature of most high-velocity gas is still poorly understood, mainly because 
the cloud distances are highly uncertain. Some of the HVCs arc clearly material stripped out 
of the Magellanic Clouds, the "Magellanic Stream" (Mathcwson, CIcary, & Murray 1974; 
Putman et al. 1998; Lockman et al. 2002), and some high- velocity absorption lines observed 
towards disk stars are most likely related to the interaction between stars/supernovae and 
the ISM (e.g., Cowie, Songaila, & York 1979; Trapero et al. 1996; Jenkins et al. 1998, 2000; 
Welty et al. 1999; Tripp et al. 2002). Apart from these clouds, direct distance constraints 
are scarce and difficult to obtain (Wakker 2001), and consequently a variety of HVC models 
remain viable. The HVCs may have a galactic origin such as the Galactic fountain (Shapiro 
& Field 1976; Bregman 1980), or they may be extragalactic (e.g., Oort 1970; Blitz et al. 
1999; Braun & Burton 1999). 

Given the difficulty of direct distance measurements, it is important to explore other 
constraints on the nature of these objects. For example, constraints on the distribution 
and size of the HVCs can be derived from studies of other galaxy groups using cither QSO 
absorption line statistics (Charlton, Churchill, & Rigby 2000) or surveys for redshifted 21 
cm emission (Zwaan & Briggs 2000; Zwaan 2001; Braun & Burton 2001; Pisano & Wilcots 

2003) . For the Milky Way HVCs, it has been suggested that the Ha emission from the clouds 
constrains their distances (e.g.. Bland- Hawthorn et al. 1998). If the gas is photoionized by 
UV flux escaping from the Galaxy, then nearby clouds should be substantially brighter in 
Ha than distant HVCs. Ha emission has in fact been detected from a variety of HVCs (e.g., 
Weiner & Williams 1996; Tufte et al. 1998, 2002; Weiner, Vogel, & Wilhams 2002), and 
detailed photoionization models place the clouds roughly 5 — 50 kpc away based on the 
observed Ha intensities (Weiner et al. 2002; Bland-Hawthorn & Maloney 2002). However, 
several observations cloud the interpretation of Ha intensities. First, the Magellanic Stream, 
which has an independently constrained distance, is much brighter in Ha than predicted by 
the photoionization models (Weiner et al. 2002; Bland-Hawthorn & Maloney 2002). Second, 
the Magellanic Stream Ha emission is spatially variable (see Figure 1 in Weiner et al. 2002), 
also contrary to the predictions of the photoionization models. Finally, O VI absorption has 
been detected in a substantial fraction of the HVCs (Sembach et al. 2000, 2002; Wakker et al. 
2002). While O VI can be produced by photoionization in large, low-density intergalactic gas 
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clouds (e.g., Tripp & Savage 2000; Tripp et al. 2001), the O VI in HVCs is almost certainly 
collisionally ionized (the cloud sizes required by photoionization models are excessive for 
HVCs, see Sembach et al. 2002). These observations suggest that collisional processes 
play an important role in the ionization of HVCs and the production of Ha emission. The 
interaction of the rapidly moving clouds with ambient magnetic fields may create instabilities 
that also play a role in the gas ionization and production of Ha emission (Konz et al. 2001). 
For these reasons. Ha distance constraints from photoionization models may need to be 
revisited. 

Ultraviolet absorption lines in the spectra of background quasars and active galactic 
nuclei (AGNs) provide another sensitive probe for the study of Galactic (as well as ex- 
tragalactic) HVCs. UV absorption lines provide detailed information on abundances and 
physical conditions in the gas, and this in turn can be compared to predictions of various 
models. For example, it is expected that Galactic fountain gas would have a substantially 
higher metallicity than infalling extragalactic gas or gas stripped from a satellite galaxy. 
Relative metal abundance patterns may also provide insights on the enrichment history of 
the HVCs; if the gas is relatively pristine, overabundances of a elements or underabundances 
of nitrogen might be observed, as seen in low-metallicity stars (McWilliam 1997) and H II 
regions (Vila-Costas & Edmunds 1993; Henry, Edmunds, & Koppen 2000). 

In this paper we use the absorption line technique to explore the nature of HVC Complex 
C, a large HVC (roughly 20° x 90°) that is at least ~5 kpc away (van Woerden et al. 1999). 
A number of Complex C abundance measurements^" have been published for a variety of 
species ranging from [N/H] = —1.94 to [Fe/H] = —0.3 (Bowen, Blades, & Pettini 1995; 
Wakker et al. 1999; Murphy et al. 2000; Richter et al. 2001; Gibson et al. 2001; CoUins et 
al. 2002). However, some of these abundances may be confused by ionization effects. Species 
such as S II or Fe II can arise in ionized as well as neutral gas leading to overestimates of 
[S/H] or [Fe/H] if no ionization correction is applied, while other species such as N I can be 
more readily ionized than H I leading to underestimation of the elemental abundance. The 
most robust species for constraining the metallicity of HVCs is O I. Like sulfur, oxygen is only 
lightly depleted by dust grains (see §6.2.1 in Moos et al. 2002, and references therein), but 
more importantly, the ionization potential of O I is nearly identical to that of H I, and O I is 
strongly locked to H I by resonant charge exchange (Field & Steigman 1971). Consequently, 
oxygen abundances based on O I and H I are relatively impervious to ionization effects 
(unless the gas is quite substantially ionized). 



^^Throughout this paper we express abundances with the usual logarithmic notation, [X/Y] = log(X/Y) 
— log(X/Y)0, and we indicate the overall metallicity with the variable Z. 
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Richter et al. (2001) have presented the first measurement of [0/H] in Complex C using 
spectra of PG1259+593 obtained with the Space Telescope Imaging Spectrograph (STIS) and 
the Far Ultraviolet Spectroscopic Explorer (FUSE); they find [0/H] = — l.OSlo'g]'. Notably, 
Richter ct al. also report that nitrogen is highly underabundant (by ~ 1 dex compared to 
oxygen), and the a elements are marginally overabundant compared to iron. These results 
suggest that Complex C is a relatively pristine extragalactic cloud plunging into the Galaxy 
for the first time. This notion has been challenged by Gibson et al. (2001), who have observed 
S 11 in several directions through Complex C and find evidence of spatial variability of the 
sulfur abundance. Sulfur abundances derived from S 11 alone are prone to ionization effects, 
as noted, but recently CoUins, ShuU, & Giroux (2003) have also reported variable metaUicity 
in Complex C based on oxygen measurements. 

Here we present new constraints on the metaUicity as well as the physical conditions, 
structure, and nature of high-velocity cloud Complex C from UV absorption lines. Our 
analysis is mainly based on high- resolution cchcllc spectroscopy of 3C 351 obtained with 
STIS, but we make use of high-resolution observations of nearby sight lines to augment the 
results with spatial information, e.g., on the transverse extent of the absorbing gas. We also 
briefly discuss intermediate-velocity gas in the 3C 351 direction. The paper is organized 
as follows: We begin with a summary of 21cm emission observations of high- velocity gas 
in the direction of 3C 351 (§2), followed by a brief summary of the STIS observations 
(§3) and the absorption line measurements (§4). We then argue that ionization corrections 
are likely to be important for these clouds, and we employ coUisional and photoionization 
models to constrain their abundances (§5). The abundances indicate that Complex C and 
the intermediate-velocity gas likely have different origins; we suggest that Complex C is an 
infalling, low-metallicity cloud while the IVC is outflowing gas from the Perseus spiral arm. 
In the discussion of our measurements, we consider the implications of the density and size 
of the absorbing gas derived from the ionization models including the conflnement of the 
gas (§6.1), and we compare 3C 351 to other nearby sight lines in the vicinity of Complex 
C (§6.2). We find evidence that the lower latitude portion of Complex C is more affected 
by interactions with the ambient medium that the higher latitude region of the cloud. We 
summarize our conclusions in §7. Throughout this paper we present spectra and velocities 
relative to the Local Standard of Rest (LSR).^^ 



^^Wc adopt the "standard" definition of the Local Standard of Rest (Delhaye 1965; Kerr & Lynden-Bell 
1986), in which the Sun is moving in the direction I = 56°, 5 ~ 23° at 19.5 km s^^. With this convention, 
the conversion to hehocentric velocity is given by fheiio — VhSR — 16.5 km s^^ for 3C 351. 
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2. High- Velocity Clouds Toward 3C 351 

The sight hne to 3C 351 (/ = 90.08°, b = +36.38°) probes an intriguing region of the 
high-velocity sky. A portion of the Hulsbosch & Wakker (1988) 21 cm emission map of 
Complex C, centered on 3C 351, is shown in Figure 1. In this figure, the color scale shows 
the LSR gas velocity while the contours indicate brightness temperature. The lowest contour 
shown corresponds to A'"(H I) = 2 x 10^^ cm~^. Complex C shows several well-defined cores; 
the 3C 351 sight hne is roughly equidistant in projection from the CIB, CD, and CK cores 
[see Wakker (2001, and references therein) for definitions of the cores and nomenclature]. CK 
has a lower velocity than most of Complex C and may be more closely associated with the 
intermediate-velocity cloud Complex K than the high- velocity Complex C cloud (see Wakker 
2001 and Haffner, Reynolds, & Tufte 2001). Figure 2 shows the intermediate-velocity 21 cm 
sky in the vicinity of 3C 351. The map in Figure 2 plots H I column densities derived 
from the 21 cm Leiden-Dwingeloo Survey (LDS; Hartmann & Burton 1997) integrated over 
—100 < i'LSR < —55 km s~^. At I > 100°, the intermediate-velocity sky is dominated by the 
"IV Arch" , but at lower latitudes it is unclear whether the emission is associated with the IV 
Arch, Complex K, or a lower- velocity region of Complex C (see Wakker 2001). We provide 
evidence below that the CK IVC and Complex C proper have different origins (§§5.4,6). 
We shall refer to the 3C 351 intermediate- velocity absorption system as C/K hereafter and 
Complex C proper (vlsr ~ —130 km s~^) as simply Complex C. 

A notable feature of Complex C is the presence of 21 cm emission at substantially higher 
velocities (ulsr < —170 km s~^) along the central ridge of the cloud superimposed on the 
main HVC emission (—155 < "Wlsr ^ —80 km s~^). The higher velocity emission is indicated 
with half-circles in Figure 1 (see also Figure 6 in Wakker 2001). The nature of the higher 
velocity component is not entirely clear, but it has a similar morphology to the lower velocity 
Complex C gas, and it is centered on Complex C. For convenience, we refer to this cloud as 
the "high- velocity ridge" (HVR) in this paper. A variety of UV absorption fines are detected 
from the high- velocity ridge in the 3C 351 spectrum as well as Complex C, and unlike C 
vs. C/K, the UV absorption lines provide evidence that the HV ridge and Complex C are 
closely related (see below). 

The H I column densities in Complex C and the high-velocity ridge are shown in Fig- 
ures 3 and 4, respectively. Figures 3-4 show 21 cm emission maps from the LDS integrated 
over —155 < i'lsr < —100 km s~^ (Complex C proper) and —205 < t;lsr < —170 km 

(high- velocity ridge). In these figures, the color scale and thin contours reflect A^(H I) 
according to the scale at the bottom. The last thin contour is drawn at the 3a limit of 
the LDS, A^(H I) = 5 x 10^^ cm"^. The thick line shows the 2 x 10^^ cm"^ contour from 
the more sensitive Hulsbosch & Wakker map in Figure 1 integrated over the same velocity 



- 7- 




i. Map of HVC Complex C 21 cm emission, from Hulsbosch & Wakker (1988), in Galactic coordinates (longitude 

increases from right to left along the x-axis). This map is centered on the sight line to 3C 351 and shows only a portion of the 
cloud (for a map of the entire HVC, see Figure 6 in Wakker 2001). Emission from high- velocity gas in the Outer Arm is also 
evident in this map closer to the plane; at Z > 90° and b < 30°, the high- velocity emission is dominated by the Outer Arm. At 
lower longitudes the distinction is less clear; the line from (/, b) = (70, 25) to (88,29) roughly delineates these structures. LSR 
velocities are indicated by color using the scale at the bottom of the figure, and the contours show brightness temperatures of 
0.05, 0.4, and 1.0 K. The half-circles show the high- velocity ridge at i;lsr ~ —200 km s~^; this higher velocity feature is clearly 
detected in absorption toward 3C351 (see text and Figures 6-8). The emission cores CIA,CIB, Cel-CelV, CK, CD, and D are 
also labeled along with other extragalactic sight lines of interest discussed in this paper. 
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Fig. 2. — Map of 21 cm emission from intermediate- velocity gas in the vicinity of 3C 351 from 
the Leiden-Dwingeloo survey including the following velocity range: —100 < wlsr ^ ~55 
km s~^. Emission in this velocity range is dominated by the "IV Arch" at / > 98°; emission 
at lower latitudes may be associated with Complex K (see §4.24 and §4.27, respectively, in 
Wakker 2001) or the Perseus Arm (see text, § 6.2.3). 



Fig. 3. — Partial map of Complex C 21 cm emission as in Figure 1, but based on data from 
the Leiden-Dwingeloo Survey (Hartmann & Burton 1997). Emission at 6 < 27° is mainly 
from the Outer Arm; higher-latitude emission is associated with Complex C. In this figure, 
colors indicate the H I column density as shown by the scale at the bottom, integrated from 
■^LSR = —155 to —100 km s~^. The lowest contour from the LDS data corresponds to the 3a 
limit of ~ 5 X 10^^ cm~^. The thick line indicates the A^(H I) = 2 x 10^^ cm~^ contour from 
the Hulsbosch & Wakker map shown in Figure 1. 
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Fig. 4. — Map of 21 cm emission from the "High- Velocity Ridge" (see text) from the Leiden- 
Dwingeloo survey. This figure is identical to Figure 3, except the emission is integrated over 
a different velocity range: —205 < fLSR < — 155 km s~^. 
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range. Much of the high-velocity emission detected closer to the plane in Figures 1-3 is 
mainly associated with the Outer Arm and the warp of the outer Milky Way (Wakker 2001, 
and references therein). 

A potentially serious source of systematic error in HVC abundance measurements is the 
large beam of the 21 cm observations usually used to estimate A^(H I). The large radio beam 
may dilute the HVC 21 cm emission, or it may show emission from gas which is in the radio 
beam but is not present along the pencil beam to the UV source. Comparisons of A^(H I) 
measurements from different radio telescopes (with different beam sizes) indicate that the 10' 
Effelsberg beam is often sufficiently small to provide a reliable H I column density (Wakker 
et al. 2001), but there are cases where even the Effelsberg beam is too large (e.g., Mrk 205, 
see Wakker 2001). Therefore it is worthwhile to review the available iV(H I) measurements 
for the 3C 351 HVCs. The HVC 21 cm emission profiles in the direction of 3C 351 observed 
with the NRAO 43m telescope (Murphy, Sembach, & Lockman 2002) and the Effelsberg 
100m telescope (Wakker et al. 2001) are shown in Figure 5. H I column densities in the 3C 
351 HVCs derived from these data are summarized in Table 1 along with iV(H I) from the 
LDS data. The uncertainties listed in Table 1 are statistical errors only. Wakker et al. (2001) 
note that the systematic uncertainties in the 3C 351 H I column densities are likely much 
larger than the statistical errors: they estimate that systematic errors lead to an uncertainty 
of ~ 1.5 X 10^^ cm~^ in A^(H I). With this uncertainty, the A^(H I) measurements in Table 1 
for Complex C are in reasonable agreement. For Complex C, we adopt A'^(H I) from the 
smallest-beam observation (Effelsberg), but with the larger systematic uncertainty reported 
by Wakker et al., i.e., N{R I) = (4.2 ± 1.5) x 10^^ cm-^. However, while both Complex C 
proper and the high-velocity ridge are apparent in the NRAO and LDS data, only Complex 
C is clearly detected in the Effelsberg observation. This suggests that the observations with 
larger beams are picking up emission which is not present along the pencil beam to 3C 351. 
Consequently, we take the NRAO Green Bank A^(H I) as an upper limit for the high- velocity 
ridge, and we derive lower limits on the HV ridge abundances below. Similarly, the H I 
column for Complex C/K is highly uncertain given the systematic uncertainty above, and 
we can only place lower limits on the metallicity of this JVC. However, these lower limits 
will turn out to be useful (§§5.3,5.4). 
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Fig. 5. — H I 21 cm emission profiles in the direction of 3C 351, recorded with (a) the NRAO 

43m telescope (Murphy, Sembach, & Lockman 2002), and (b) the Effelsberg 100m telescope 
(Wakker et al. 2001). Brightness temperature is plotted vs. LSR velocity. Complex C 
proper and the high-velocity ridge are marked along with the intermediate-velocity cloud 
Complex C/K. HI column densities derived from these data are listed in Table 1. 
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Table 1. Intermediate- Velocity and High- Velocity Cloud H I Column Densities toward 3C 

351 





Efrclril)(>rg 


lOOm (10' bourn)" 


GrcQii Buuk 43in (21' 1k>;uii)'' 


D\viiig(>loo 


25m (35' ]>Qiimy 


High- Velocity Cloud 




7V(H I)'^ 


VhSR 


N{R l)'^ 


VhSR 






(km s-i) 


(IQis cm-2) 


(km s-i) 


(IQis cm-2) 


(km s-i) 


(10^8 cm-2) 


High- velocity ridge .... 






-182.2 


4.19 


-191 


4.5 ±0.9 


Complex C (proper) . . . 


-129 


4.2 ±0.3 


-128.2 


6.01 


-126 


4.5 ±0.6 


Complex C/K 


-89 


1.3 ±0.2 


-83.0 


4.08 


-69 


4.0 ±0.5 



^From Wakkcr ct al. (2001). 
•"From Lockman ot al. (2002). 

'^From Wakker et al. (2001), derived from the Lcidcn-Dwingeloo Survey (Hartmann & Burton 1997). 

'^Listed column density uncertainties include only statistical error, and uncertainties are not reported by Lockman 
et al. (2002) for the Green Bank 43m measurements. Wakker et al. (2001) discuss several sources of systematic error, 
and they estimate that these lead to an uncertainty of ~ 1.5 x 10^^ cm~^ in A''(H I). 
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Fig. 6. — Portion of the FUV spectrum of 3C 351 obtained with STIS in the E140M echelle 
mode (FWHM ^ 7 km s~^) plotted versus observed wavelength. For display purposes, the 
spectrum has been binned (2 pixels into 1) in this figure (all other figures, as well as the 
line measurements, make use of the full resolution, unbinned spectrum). This region of 
the spectrum shows the absorption profiles of the O I A1302.2 and Si II A1304.4 lines due 
to the ISM of the Milky Way as well as several extragalactic absorption lines. Four main 
components are readily apparent in the Galactic O I and Si II profiles; these are indicated 
with tick marks immediately above the spectrum. We are mainly interested in the high- 
velocity absorption lines associated with Complex C and the "high-velocity ridge" (HVR, 
see §2), which are marked with longer tick marks. Absorption arising in the intermediate- 
velocity cloud Complex C/K and local gas is also identified with short ticks. The thin 
solid line indicates the fitted continuum with the predicted Lj6 lines due to the associated 
absorption systems at ;2abs = 0.3709 and 0.3719 superimposed (see Yuan et al. 2002). The 
other extragalactic lines in this wavelength range are identified as Lya at Zabs = 0.0699 and 
0.0715 (see text for further details). 
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Fig. 7. — Effelsberg 21 cm emission profile (top left panel) and selected UV absorption line 
profiles of Milky Way lines detected toward 3C 351, plotted versus Local Standard-of-Rest 
(LSR) velocity. Brightness temperature is plotted in the top left panel, and observed flux is 
shown in the other panels. Dashed lines show the fitted continua and dotted lines mark the 
flux zero level. Absorption lines due to Complex C and the "High- Velocity Ridge" (HVR, 
see text) are indicated with tick marks. Unrelated absorption features are also marked (x). 
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3. STIS Echelle Spectroscopy 

We now turn to the UV absorption line data. The STIS observations of 3C 351 are fully 
described in §2 of Yuan et al. (2002). Briefly, the QSO was observed with the E140M echelle 
mode with the 0'.'2 x 0'.'06 slit; this mode provides 7 km s~^ resolution (FWHM) and covers 
the 1150 - 1710 A range with only a few gaps between orders at the longest wavelengths 
(see Woodgate et al. 1998 and Kimble et al. 1998 for further details). The data were 
reduced following standard procedures with software developed by the STIS Investigation 
Definition Team. Scattered light was removed using the procedure developed by the STIS 
Team (Bowers et al. 1998), which corrects for echelle scatter as well as other sources of 
scattered hght. 

A sample of the final spectrum is shown in Figure 6. This somewhat complicated portion 
of the spectrum shows the absorption profiles of the important Milky Way O I A1302.2 and 
Si II A1304.4 profiles as well as several extragalactic absorption lines. Absorption profiles of 
several species of interest are plotted versus LSR velocity in Figure 7 along with the H I 21 
cm emission observed in the direction of 3C 351 by Wakker et al. (2001) with the 100 m 
Effelsberg telescope. The UV absorption lines show four distinct components. The two 
highest- velocity components at wlsr = —128 and —190 km s~^ are close to the observed 
21 cm velocities of Complex C and the high- velocity ridge. Similarly, the UV fines at vlsr 
«i —82 km s~^ are close to the velocity of IVC C/K. Analysis of the Complex C/K absorption 
lines is complicated by substantial blending with adjacent components. Nevertheless, the 
data indicate that IVC C/K has a higher metalficity than HVC Complex C (§5.4). 

In this paper, we are primarily interested in Complex C and the high-velocity ridge. 
However, absorption from other Galactic and extragalactic clouds can introduce systematic 
errors that must be considered. For example, 3C 351 has a dramatic "associated" absorption 
system (i.e., z^hs ~ ^Qso) that affects the spectrum near the N V, O VI, and Lyman series 
lines at tfie redsfiift of tfie QSO (Yuan et al. 2002, and references tfierein). Tfie associated 
Lj6 lines at z^^y^g = 0.3709 and 0.3719 fall near tfie Galactic O I A1302.2 profile, wfiicfi is 
a crucial profile for our analysis. We can predict tfie strengtfi of tfiese contaminating Lj6 
fines based on fits to the other Lyman series fines in the STIS spectrum from Yuan et al. 
(2002). The thin solid fine in Figure 6 shows the predicted LyS lines superimposed on the 
fitted continuum. Fortunately, as can be seen from Figure 6, these associated LyS lines have 
no impact on the high- velocity O I lines of primary interest. Tfie nearby Lya lines at Zabs 
= 0.0715 and particularly 0.0699 could be a more serious source of contamination tfian tfie 
associated LyS lines. Tfiese Lya lines appear to be mostly separated from tfie Milky Way 
O I and Si II profiles, but since Lya lines are clustered on velocity scales of a few fiundred 
km s"^ (Penton et al. 2000), it remains possible tfiat unseen Lya lines clustered around tfie 
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strong Zahs — 0.0699 and 0.0715 absorbers affect the Galactic O I and/or Si II profiles in 
Figure 6. However, a detailed comparison of the Milky Way profiles in the following section 
indicates that there is little or no contamination present in the O I and Si II lines. 

The Galactic Si II Al 193.3 transition is badly blended with C III A977.0 absorption from 
a Lyman hmit system at Zghs — 0.2210, so we do not use this transition for the HVC analysis. 
Also, the Milky Way Si II A1190.4 and S III A1190.2 transitions are highly blended, and we 
do not trust information from either of these lines. Finally, we note that the Galactic Si II 
A1260.4 and C II A1334.5 transitions are severely saturated and consequently do not provide 
useful constraints. 



4. Absorption Line Measurements 

To measure the line column densities and to assess the impact of unresolved saturation 
on the absorption lines of interest, we mainly rely on the apparent column density technique 
(Savage & Sembach 1991), but some supporting measurements are also derived from Voigt 
profile fitting, using the software of Fitzpatrick & Spitzcr (1997) with the line spread func- 
tions from the STIS Instrument Handbook (Leitherer ct al. 2001). The apparent column 
density is constructed using the following expression, 

N,{v) = {m,c/ne^){fX)-'T,{v) = 3.768 x 10^^(/A)-Hn[/e(^;)//(^;)]. (1) 

This provides the instrumentally broadened column density per unit velocity [in atoms cm~^ 
(km s"^)""-*^], where / is the transition oscillator strength, A is the transition wavelength 
(the numerical coefficient is for A in A), I{v) is the observed line intensity and Idv) is the 
estimated continuum intensity at velocity v. By comparing two or more resonance lines of a 
given species with adequately different /A values, the profiles can be checked for unresolved 
saturation and a correction can, in some cases, be applied (Jenkins 1996). If the Ng,{v) profiles 
of the resonance lines are in good agreement, then the lines are not affected by unresolved 
saturation, and the profiles can be integrated to obtain the total column, A^tot = / Ni^{v)dv. 

High-velocity absorption line equivalent widths and apparent column densities, inte- 
grated across the velocity ranges that delimit Complex C proper and the high- velocity ridge, 
are summarized in Table 2. These measurements were made using the methods of Sembach 
& Savage (1992) and include continuum placement uncertainty and a 2% flux zero point un- 
certainty in addition to the statistical uncertainties. Golumn densites of well-detected species 
determined from Voigt proflle fltting (Apf) are also listed. Equivalent widths and column 
densities for the intermediate-velocity G/K cloud are provided in Table 3. The abundance 
measurements presented in the last two columns of Tables 2-3 are discussed in §5. 



Table 2. High- Velocity Absorption Lines toward 3C 351: Complex C and the 

High- Velocity Ridge 



Species 




r 


Cloud 


<v >^ 




log iVa^ 


log iVpfd 


[XVH 1]'^ [X/H]f 




(A) 






(km s""'^) 


(mA) 









1 1302.17 0.0519 Complex C 

HV Ridge 

N 1 1199.55 0.130 Complex C 

HV Ridge 

Si Ift .... 1304.37 0.0917 Complex C 

HV Ridge 
1526.71 0.132 Complex C 
HV Ridge 

AlH 1670.79 1.83 Complex C 

HV Ridge 

SII 1259.52 0.0166 Complex C 

HV Ridge 

FeH 1608.45 0.0580 Complex C 

HV Ridge 

Si HI.... 1206.50 1.67 Complex C 

HV Ridge 

Si IV 1393.76 0.514 Complex C 

HV Ridge 

CIV 1548.20 0.191 Complex C 

HV Ridge 



-128 ± 1 
-190 ±2 

-130 ±3 



-128 ± 1 
-183 ±5 
-128 ± 1 
-184 ±3 

-127± 1 
-182 ±3 

-124 ±6 

-125 ±2 



-127 ±2e 
-189 ±3 



115.1±4.9 
81.8±8.1 

(24.6±10.8) 
(3.7±14.6) 

63.3±4.6 
30.3±7.1 
86.1±6.5 
68.8±10.5 

105.0±6.9 
70.1±11.6 

(12.0±5.6) 
(1.1±7.9) 
63.1±8.7 

(12.4±13.8) 

115.2±7.6' 
196.0±12.5 



14.45±0.04 
14.12±0.05 

<13.44S 
<13.55« 

13.81±0.05 

13 42+0 08 

13.72±0.05 

1 Q 4Q+0 06 



12.57±0.04 
1 2 sn+0 06 

iz.ou_o.o8 



< 14.08 

< 14.1S 



14.63 



-1-0.10 
-0.14 



-0.73 



-1-0.23 
-0.21 



13.78 



-0.06 
-0.08 



< 13.6S 



14.17±0.04 > -1.24 



< -1.11 



13.77±0.04 
13.50±0.05 
13.77±0.04 
13.50±0.05 

12.56±0.05 
12.33±0.05 



-0.40t°;?° 

> --0.70 
-0 40+0-20 

> -0.70 



-0 52+0-21 
> -0.79 



< 0.05 



^•'0-0.21 

> -1.26 



< -1.37 



n 70-1-0.20 
^•'"J-0.14 

> -1.26 

r> 70-1-0.20 

^•'•'-O.M 

> -1.26 



-0. 
> -1.36 



0.21 

0.15 



< -0.31 



13. 



f0.14 
-0.21 



-0.29tg-j^ -0.5410-22 



12.91±0.06' ••• 

> 13.07j 13.40±0.10j 



-118 ±13 (7.1±6.2) < 12.4s 
-188 ±4 45.1±8.4 



12.771^^ 12.79±0.07 



-120 ±8 (14.6±8.9) < 13.08 
-191 ±3 100.2±11.6 13.50±0.06 13.54 



-fO.16 
-0.26 



Table 2 — Continued 



Species Aq^ f Cloud < v >^ Wx" log N^" log Npi'^ [X'/H I]'' [X/H]f 
(A) (km s-i) (mA) 



"^Rest-frame vacuum wavelength and oscillator strength from Morton (2002) or Morton (1991). 
'^Profile- weighted mean velocity of the line, < v >= / v[l — I{v) / Ic{v)]dv / /[I — I{v) / Ic{v)\dv. 

'^Integrated equivalent width and apparent column density (in cm~^). For the High- Velocity 
Ridge, all quantities are integrated from i^lsr = —230 to —155 km s~^, and the Complex C 
absorption lines are integrated over —150 < v\,sk < —110. Equivalent widths in parentheses have 
less than 3c7 significance, and we do not consider these lines to be reliably detected. 

''Column density (in cm~^) from Voigt profile fitting. 

'^Implied logarithmic abundance if ionization corrections are neglected, i.e., [X'/H I] = log 
7V(X*)/7V(H I) - log (X/H)q. 

Logarithmic abundance obtained by applying the ionization correction from the best-fitting models presented in 
the text. In this table we have used the collisional ionization equilibrium corrections, but very similar results are 

obtained from the CLOUDY model (sec§§ 5.2.1,5.2.2,5.3). Error bars include column density uncertainties and solar 
reference abundance uncertainties but do not reflect uncertainties in the ionization correction. 

^4(7 upper limit assuming the linear curve-of-growth applies. 

'^The Voigt profile fitting column is from a simultaneous fit to the 1304.37 and 1526.71 A lines. 

Similarly wc combined the information from these two transitions for the directly integrated column. 
For the final integrated Si II column density, wc adopt the following weighted averages of the 1304.37 
and 1526.71 A measurements. Complex C: log N{S\ II) = 13.76±0.03; high-velocity ridge: log 
7V(Si II) = 13.46±0.06. 

'Formal error bars may underestimate the true uncertainty due to strong blending with lower 
velocity gas (see Figure 7). 

j Strongly saturated absorption. 
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Table 3. Intermediate- Velocity Absorption Lines toward 3C 351: Complex C/K^ 



Species 


Ao 


< V > 




log 


log A^pf 


[XVH I]b 






(A) 


(km s"^) 


(mA) 












I 


1302.17 


-82 ± 1 


180.6±5.4<^ 


> 14.66'' 




> -0.50 


> 


-0.52 


N I 


1199.55 


-83± 10 


66.5±13.5 






> -0.91 


> 


-0.69 


Si II 


1304.37 


-81 ±1 


131.3±4.9 


> 14.23** 




> 0.20 


> 


-0.52 




1526.71 


-81 ± 1 


181.1±7.0 


> 14.09"= 


14 oc+0.13f 


> 0.20 


> 


-0.52 


Al II 


1670.79 


-82 ± 1 


191.1±7.2d 


12.88±0.04'i'S 


~ 13.1^ 


> -0.20 


> 


-0.88 


S II 


1259.52 


-86 ±4 


26.1±6.2 


i4.io;°:i° 




> 0.16 


> 


-0.46 


Fe II 


1608.45 


-84 ±2 


79.9±10.2 


13.91±0.07 


13.86±0.10 


> -0.15 


> 


-0.64 


Si IV 


1393.76 


-83 ±2 


68.9±6.3 


12.98±0.05 


13.03±0.06 








C IV 


1548.20 


-82 ± 1 


110.6±9.0'" 


13.58±0.05'^ 


13.29±0.07'' 









"^See Table 2 for definitions of the quantities in this table. The equivalent widths and apparent 
column densities are integrated over the velocity range of the intermediate- velocity C/K cloud, from 
^^LSR = —110 to —59 km s~^. 

^Implied logarithmic abundance if ionization corrections are neglected, i.e., [X*/H I] = log 

N{X^)/N{R I) - log (X/H)0. 

"^Logarithmic abundance obtained by applying the ionization correction from the best-fitting 
CLOUDY model presented in the text (§ 5.4). Error bars include column density uncertainties 
and solar reference abundance uncertainties but do not reflect uncertainties in the ionization cor- 
rection. 

"^Line is strongly blended with adjacent absorption features. 
"Saturated absorption line. 

^Simultaneous fit to the Si II A1304.37 and A1526.71 lines. 

^Apparent column density may underestimate the true column due to unresolved saturation. 
^Poorly constrained result due to strong blending. 
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As noted in §1, the most important species for metallicity measurements is O I. The only 
O I transition available in the direction of 3C 351 is the 1302.2 A line.^^ In the high-velocity 
ridge, N{0 I) appears to be well-constrained from the 1302.2 A profile. In Complex C, on 
the other hand, the 1302.2 A line does not go to zero flux in the core, but it is quite strong 
(see Figs. 6-7). Consequently, it is important to consider whether the Complex C O I column 
might be underestimated due to unresolved saturation. We believe that the O I line is not 
seriously saturated based on the following evidence. We first note that the Si II hues show 
no indications of unresolved saturation. Figure 8a compares the Na,{v) profiles of the 1304.4 
and 1526.7 A lines of Galactic Si II, transitions that differ in /A by 0.23 dex. These Si II 
Ng_{v) profiles are in good agreement, ■'^^ which indicates that the hues are not significantly 
saturated. If O I is affected by unresolved saturation, then the O I Ng^{v) profile should differ 
from the Si II Ni^{v) profiles due to the underestimation of A^(0 I) in the saturated pixels. 
Figure 8b compares the O I apparent column density to a weighted composite Si II Na{v) 
profile constructed from the 1304.4 and 1526.7 profiles following the procedure of Jenkins 
& Peimbert (1997), with the O I profile scaled down by a factor of five. From Figure 8b, 
we see that the O I and Si II profiles have very similar shapes. This suggests that there 
is little unresolved saturation within the O I profile. We notice from Table 2 that Voigt 
profile fitting gives a somewhat higher A'"(0 I) than direct Na{v) integration, although the 
O I columns from the two methods agree within the la uncertainties. This could refiect a 
small amount of saturation. To be conservative, we adopt the higher A'^(0 I) value for our 
subsequent analysis. However, we will consider the lower O I column in our analysis, and we 
will find that our conclusions are essentially the same. 

The comparison of the O I and Si II Ng^{v) profiles in Figure 8b also indicates that there 
is little contamination of the O I lines from the nearby Lja absorbers discussed in §3 and 
shown in Figure 6. Significant optical depth from such contamination would cause the O I 
and Si II profiles to have different shapes at velocities where the contamination occurs. The 
O I and Si II profiles have similar velocity structure, and moreover the O I/Si II column 
density ratio is very similar in Complex C and the high-velocity ridge. The Al II and Si II 
Ng_{v) profiles are also remarkably alike, as shown in Figure 8d. The Fe II Na,{v) profile is 
compared to the composite Si II profile in Figure 8c. Fe II is clearly detected in Complex C, 



^^In principle, weaker O I transitions can be observed with FUSE at 915 < A < 1050 A. Unfortunately, 

there is not sufficient flux from the QSO in this wavelength range for absorption line measurements because 
the spectrum is strongly attenuated below 1115 A by an optically thick Lyman limit absorber at Zabs = 
0.2210 (see Mallouris et al., in preparation). 

^^Unresolved saturation would be manifested by a higher apparent column in the weaker transition (1304.4) 
compared to the stronger transition. The highest outlying point in the Na,{v) profiles in Figure 8a is from 
the 1526.7 A transition, the opposite of the expected effect. 
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but the substantially greater noise limits comparisons of C and the HV ridge for this species. 
The similarity of the O I, Si II, and Al II profiles in Figure 8 would seem to suggest that the 
physical conditions and relative abundances in Complex C and the high-velocity ridge are 
roughly the same. However, the behavior of the higher ionization stages (Si IV and C IV) is 
entirely different in Complex C and the high-velocity ridge, as shown in in Figure 8e-f. Very 
little high-ion absorption is apparent in the velocity range of Complex C, but Si IV and C IV 
are clearly detected in the high-velocity ridge, with component structure similar to that of 
the lower ionization stages. The Si IV and C IV profiles therefore indicate that there is a 
important difference in the physical conditions of Complex C and the HV ridge. We discuss 
this issue further in the following section. 
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Fig. 8. — Apparent column density as a function of LSR velocity {Ng,{v), see text § 4) of 
selected absorption lines detected in the high- velocity clouds in the direction of 3C 351. (a) 
Na,{v) profiles of the Si II transitions at 1304.4 A (filled circles) and 1526.7 A (open circles). 
The good agreement of these transitions indicates that the fines are not strongly affected by 
unresolved saturation, (b) O I A1302.2 profile scaled by 0.2 (open squares) vs. a weighted 
composite Si II profile (histogram), (c) Fe II A1608.5 profile (filled circles) vs. the composite 
Si II profile (histogram), (d) Al II A1670.8 profiles scaled by 17.0 (open squares) vs. the 
composite Si II profile (histogram), (c) Si IV A1393.8 profile scaled by 5.0 (filled circles) 
vs. the composite Si II profile (histogram), (f) C IV A1548.2 profile (open squares) vs. the 
composite Si II profile (histogram) . 
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Fig. 8. — continued 
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Fig. 8. — continued 
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5. Abundances and Ionization 

We next consider the implications of the measurements presented above. We first argue 
that ionization corrections are important in the 3C 351 HVCs (§5.1). We then investigate 
the absolute and relative abundances in Complex C and the high- velocity ridge with the aid 
of ionization models (§§5.2—5.4). 

5.1. Preliminciry Renicirks 

In general, the logarithmic abundance of species X with respect to species Y is 

where N(X.^) and /(X') are the column density and ion fraction of the i*^ ionization stage 
of species X (and likewise for Y), and (X/Y)© is the solar reference abundance. In many 
situations, the ionization correction /(Y*)//(X*) can be neglected, e.g., if /(Y*) = /(X*) = 1. 
However, when we state that the "ionization correction is decreasing" , this does not indicate 
that the ionization correction is becoming negligible; this only indicates that /(Y*)//(X*) 
is decreasing. In this paper, to distinguish between measurements that have been corrected 
for ionization and those which have not, we adopt the following notation: uncorrected abun- 
dances are indicated by the observed ion, e.g., [Si II/H I], while measurements which have 
had an ionization correction applied are indicated without reference to a particular ion, 
e.g., [Si/H]. The penultimate columns of Tables 2-4 list the uncorrected abundances im- 
plied by our column density measurements. The last columns of these tables summarize 
the ionization-corrected abundances, using the ionization models described in the sections 
below. For the high- velocity ridge and IVC C/K, we list lower limits on the abundances 
because we only have upper limits on A^(H I) for these components (see § 2). Ionization 
corrections can be a large source of uncertainty; we present examples below. The A^(H I) 
measurements are often the other main source of uncertainty. We include our best estimates 
of the A'^(H I) uncertainty in our abundance error bars. However, systematic errors in A^(H I) 
due to radio beam effects can be much larger than statistical errors, and these uncertainties 
can be difficult to assess. 



^^We adopt the solar abundances reported by Holweger (2002) for the most abundant elements. The oxygen 

abundance from Holweger is in excellent agreement with the independent solar oxygen measurement reported 
by AUende Prieto, Lambert, & Asplund (2001), and these solar abundances are close to the interstellar oxygen 
measurements in the vicinity of the Sun (Sofia & Meyer 2001). Solar abundances of sulfur and aluminum 
are taken from Grevesse, Noels, & Sauval (1996). 
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If we neglect the ionization correction for the Complex C absorption lines toward 3C 
351, we obtain highly discrepant results from neutral-gas tracers compared to species that 
can persist in ionized gas, i.e., the different ions imply significantly different metallicities. For 
example, we obtain [O I/H I] = -0.7 and [N I/H I] < -1.1 vs. [Si II/H I] = -0.4, [Fe II/H I] 
= —0.3, and [Al II/H I] = —0.5. If correct, these abundances would indicate a very unusual 
relative abundance pattern. It is more likely that these discrepancies indicate that ionization 
corrections are important for the 3C 351 sight line. These abundance discrepancies are larger 
than the uncertainties in the respective measurements, as shown in Figure 9. In this figure 
we plot relative abundance measurements or limits, with respect to Si II, for Complex C, the 
high- velocity ridge, and IVC C/K derived directly from the O I, N I, Fe II, Al II, and S II 
lines, i.e., without ionization corrections. We plot abundances relative to Si II because we do 
not have a robust A^(H I) measurement in the high- velocity ridge or IVC C/K, and because 
Si II is the best-constrained species in Complex C. The abundances in Figure 9 are arranged 
according to ionization potential. Again, we see evidence of ionization effects in all three of 
these clouds: the neutral species O I and N I are significantly underabundant with respect 
to Si II while ionized-gas tracers (Fe II and Al II) are much closer to the solar pattern. 

5.2. Complex C Proper 

5.2.1. Collisional Ionization 

Can we reconcile the various abundances by applying ionization corrections? Since we 
have reasons to believe that collisional processes may be important (see §1), and in order to 
make our development more clear, we begin with a simplified picture where coUisional ioniza- 
tion equilibrium applies, with no additional photoionization from an energetic radiation field 
(we will add photoionization in the next section). Figure 10 shows relevant column density 
ratios, as a function of temperature, derived from the equilibrium collisional ionization calcu- 
lations of Sutherland & Dopita (1993) with solar reference abundances from Holwcgcr (2001) 
and Grevesse et al. (1996).^^ The observed ratios in Complex C arc ovcrplottcd for com- 
parison. We see that over the temperature range in which the O l/Si II ratio is reproduced 
(within Icr), the model does not exactly match the other observed ratios; departures from 
solar relative abundances are required. Nitrogen, for example, must be made underabundant 



^^Most elements in the Sutherland & Dopita (1993) calculation begin in the dominant ionization stage in 

H I regions at the lowest temperature considered. However, iron is an exception; the Fe II ion fraction is 
significant at the lowest T. We have assumed that Fe starts from Fe II in the ionization ladder because any 
realistic radiation field will fully ionize the Fe I. 
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Fig. 9. — Logarithmic abundances, with respect to silicon, of O I, N I, Fe II, Al II, and S II 
in Complex C (filled circles), the high- velocity ridge (open squares), and the intermediate- 
velocity cloud C/K (open triangles) in the direction of 3C 351. The abundances in the 
outer arm component toward H1821+643 are also shown (open diamonds). The species are 
plotted in order of increasing ionization potential (shown in parentheses in eV below each 
element; the IP of Si II is 16.35 eV). No ionization corrections have been applied. The 
underabundances of O I and N I with respect to Si II suggest that ionization corrections are 
important. 
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with respect to silicon by at least 0.5 — 0.7 dex (an underabundance moves the model curves 
down in Figure 10). This may not be surprising since there is other evidence that nitrogen is 
underabundant in Complex C (Richter et al. 2001; Collins et al. 2003), and nucleosynthesis 
in low-metallicity gas can result in N underabundances (e.g., Henry et al. 2000). Similarly, 
Figure 10 indicates that aluminum must be underabundant if the gas is collisionally ionized 
and in equilibrium. Like nitrogen, Al could be underabundant, although much less so, due 
to nucleosynthesis effects since it is an odd-Z element, as observed in low-metallicity Galac- 
tic field stars (e.g., Lauroesch et al. 1996). An Al underabundance in the gas phase could 
alternatively be due to depletion by dust since aluminum is highly refractory (e.g., Jenkins 
1987). However, this would be inconsistent with the iron abundance: Fe II should also be 
depleted in this case, but Figure 10 seems to require a slight Fe overabundance, by ~0.2 dex, 
which leaves little room for dust depletion. 

The only puzzhng imphcation of Figure 10 is, in fact, the Fe overabundance. Nitrogen 
and Al underabundances would imply, at face value, that the gas is relatively pristine. 
However, in this case Fe should be underabundant as well since iron is thought to be mainly 
synthesized on longer timescales in Type la supernovae. We note the Fe II profile is the 
noisiest measurement among the detected lines (compare Figure 8c to the other panels in 
Fig. 8), and the observed Fe Il/Si II ratio is only 2a above the model prediction at log T ~ 4.2. 
The shght Fe overabundance may simply be a result of noise. In fact, we shall see below that 
when uncertainties in the reference abundances as well as the column density measurements 
are taken into account and ionization corrections are applied, the final abundances agree 
within 1(7, except nitrogen. However, it is clear that iron is not highly underabundant due 
to dust depletion or nucleosynthesis effects. We also note that Murphy et al. (2000) find a 
relatively high Fe II/H I ratio for Complex C in the direction of Mrk 876. For Figure 10, we 
assumed solar relative abundances. If we were to adopt overabundances of the «— elements 
(O I, Si II, and S II) by ~0.3 dex, on the grounds that such patterns arc observed in low- 
metallicity stars (McWilliam 1997, and references therein), then we would have a substantial 
discrepancy between the observed and model Fe II columns. 

For the calculation of abundances using equation 2, we present in Figure 11 the ioniza- 
tion corrections obtained from the coUisional ionization equilibrium calculations of Suther- 
land & Dopita (1993). The O l/Si II ratio suggests that log T ^ 4.20; taking the ion fractions 
from Figure 11 at this temperature, with the column densities from Table 2 and A'^(H I) = 
(4.2 ± 1.5) X 10^^ cm~^, we obtain the following abundances for Complex C proper: 

[0/H]c = -0.76l°i?, (3) 

[N/H]c < -1.37, (4) 
[Si/H]c = -0.731°:?°, (5) 
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Fig. 10. — Column density ratios predicted in collisional ionization equilibrium vs. gas 
temperture, based on the calculations of Sutherland & Dopita (1993) with solar relative 
abundances from Holweger (2001) or Grevesse et al. (1996). Ratios of various species to 
iV(Si II) are shown with small symbols according to the key on the right side. Observed 
ratios are plotted at log T = 4.2 with larger symbols and la error bars. Arrows indicate 4a 
upper limits. 
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log [T (K)] 

Fig. 11. — Ionization corrections for collisionally ionized gas in equilibrium from the calcu- 
lations of Sutherland & Dopita (1993) vs. log T. As in Figure 10, the species corresponding 
to each curve is indicated by the key at the right. The H I ion fraction (/ = H I/Htotai) is 
shown with a heavy black line using the scale on the right axis. The Al II and S II curves are 
very similar because over this temperature range, /(Al II) ^ /(S II) ~ 1, so the ionization 
correction is nearly equal to /(HI). 
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[Fe/H]c = -0.541°:??, (6) 

[Al/H]c = -O.SSt'oil, (7) 

and 

[S/H]c < -0.31. (8) 

The error bars in these abundances include the uncertainties in the column densities and the 
solar reference abundances but do not include the uncertainties in the ionization correction. 
We note that application of the ionization correction removes the abundance discrepancies 
discussed above. After the ionization correction at log T = 4-^ has been applied, all of the 
abundances agree within the la uncertainties, with the notable exception of nitrogen. The 
implied metallicity is Z — 0.1 — 0.3 Zq. The nitrogen underabundance suggests that the gas 
has undergone relatively few cycles of metal enrichment from stars. The magnitude of the 
nitrogen discrepancy decreases if we use the lower value for N{0 I) from N.^{v) integration 
(see Table 2), and the best estimate of [0/H] decreases as well, by ~ 0.2 dex. However, since 
we only have an upper limit on A^(N I), we would still be faced with a significant nitrogen 
deficit. It would be valuable to obtain additional observations of 3C 351 to more securely 
constrain the nitrogen situation (see additional discussion in § 6.2.1), and to measure the Fe 
abundance with less noise. 

The combination of low overall metallicity and a nitrogen underabundance argue that 
Complex C is not ejecta produced in a Galactic fountain, which would produce gas with a 
higher metaUicity. One might suppose that the HVC could still be part of a Galactic foutain 
if during the foutain cycle the gas is substantially diluted with low-metallicity gas or if the 
fountain flow started at large Galactocentric radii, as suggested by Gibson (2002). However, 
both of these hypotheses have trouble explaining the nitrogen underabundance, because N 
is not underabundant in the disk ISM (e.g., Meyer, Cardelli, & Sofia 1997), even at larger 
Galactocentric radii (Afilerbach, Churchwell, & Werner 1997). It seems more likely that 
Complex C has an extragalactic origin. It could be stripped gas from a satellite galaxy, or 
it could have a more distant origin. 



5.2.2. Photoionization 

Although there is strong evidence that collisional processes are important in HVCs (see 
§1), it is worthwhile to investigate how additional photoionization could alter the observed 
column densities. In addition, it remains possible that the Ha emission and high-ion absorp- 
tion arise from a coUisionally ionized phase (e.g., an interface on the surface of the HVC) 
while the low ionization lines originate inside the cloud where the gas is mainly photoionized. 
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For this purpose, we have constructed photoionization models using CLOUDY (v94.0; 
Ferland et al. 1998). The character of the radiation field to which an HVC is exposed is 
not entirely clear. This depends on the location of the HVC and the fraction of the ionizing 
photons that escape from the Calactic disk (Weiner et al. 2002; Bland-Hawthorn & Maloney 
2002). However, the extragalactic UV background from quasars and active galactic nuclei 
provides a floor; additional photons from local sources only add to this background. Conse- 
quently, we begin with the UV background from QSOs-|-AGNs only. We model the Complex 
C absorber as a plane-parallel, constant-density slab exposed to the QSO background at 
z ^ from Haardt & Madau (1996) with = 1 x 10"^^ ergs s'^ cm.-^ Hz"^ sr"^ at 1 
Rydberg (see Weymann et al. 2001 and references therein for observational constraints on 
the UV background intensity). The absorber thickness is adjusted to reproduce the observed 
A^(H I), and the metallicity and ionization parameter f/ (= H ionizing photon density/total 
H number density) are varied to match the metal column densities. It is important to note 
that the H I column is high enough so that self-shielding effects are important. These models 
should not be simply scaled for use on higher (or lower) A^(H I) absorption systems. 

Of course, these models are necessarily simplifled compared to a real absorption system. 
However, a more detailed treatment of radiation transfer, geometry, and multiphase effects 
in the model presented by Kepner et al. (1999) in the end predicts very similar column 
densities to an analogous CLOUDY model. While this is only one test case, and more 
detailed modeling is highly warranted, the good agreement with CLOUDY is encouraging. 

Figure 12 shows the metal column densities predicted by the CLOUDY model as a 
function of U, and Figure 13 shows the corresponding ionization corrections from the same 
model over the same ionization parameter range. The best fit is obtained with log U ^ —4.45, 
as shown in Figure 12. At this value of U, the derived abundances are nearly identical to 
those derived above from the coUisional ionization model. The only substantial difference 
between the photoionized and coUisionally ionized models is for nitrogen. In the coUisionally 
ionized model, all of the ionization correction factors of interest decrease as T increases 
and the gas becomes more highly ionized (see Figure 11). In the photoionized model, on 
the other hand, the N ionization correction increases while the other ionization corrections 
decrease as the gas becomes more highly ionized (see Figure 13); this is due to the relatively 
high photoionization cross section of N I. This reduces the nitrogen underabundance implied 
by the measurements somewhat compared to the coUisional model. However, the best- 
fitting photoionized model still suggests an N underabundance: we find [N/H]c < —1.01 
vs. [0/H]c = — 0.75lo;29- If take the lower A'"(0 I) from direct integration (see Table 2 
and §4), then the best fit is obtained with a somewhat higher ionization parameter, log 
U ~ —3.95 (log rzH ~ —2.55). As in the coUisionally ionized model, this requires a lower 
[0/H] by ~ 0.2 dex, and the nitrogen deficiency is reduced. The iron problem is also present 
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Fig. 12. — Model of gas photoionized by the extragalactic UV background at 2; ~ with 
log A^(H I) = 18.62, Z = O.I8Z0, and relative abundances following the solar pattern from 
Holweger (2002). Model column densities are plotted with small symbols, shown in the key 
at the right, as a function of the ionization parameter U (lower axis) and the particle density 
(middle upper axis). The mean gas temperature is shown on the uppermost axis; note that 
the temperature scale is not linear. Observed column densities are indicated with larger 
symbols with la error bars. Points with arrows are Aa upper limits. 
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Fig. 13. — Ionization corrections for the photoionization model shown in Figure 12. As in 
Figure 11, the curves are identified at right, and the thick black line indicates /(H I) using 
the right axis. 
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in the photoionized calculation: if we make a— elements overabundant in the model, then we 
find that the observed A^(Fe II) is substantially greater than the predicted column density. 
This would be unexpected in low-metallicity gas, especially if nitrogen is under abundant. 

The similarity of the photoionization and coUisional ionization results is not necessarily 
surprising. CLOUDY includes coUisional processes. The gas temperature in the photoion- 
ization model is governed by the balance of photohcating and various sources of cooling, and 
with the right combination of temperature and density, the coUisional processes may domi- 
nate. The different behavior of N in the photoionized vs. collisionally ionized models is also 
not surprising. Nitrogen has a relatively large photoionization cross section, and the nitrogen 
+ hydrogen charge exchange reaction is much weaker than that of oxygen. Consequently, 
N I is more readily ]?/iotoionized than most species, and N I is a sensitive indicator of par- 
tially ionized gas (Sofia & Jenkins 1998). We experimented with other radiation fields in the 
photoionization model with various amounts of stellar flux added to the QSO background, 
and we found that the low-ion results did not change dramatically. 



5.3. High- Velocity Ridge 

The low-ion column density ratios in Complex C proper and the high- velocity ridge are 
strikingly ahke (see Figures 8-9 and Table 2). This suggests that the physical conditions 
and abundances in Complex C and the HVR are quite similar. As discussed in §2, the 
H I column in the high-velocity ridge is not securely measured, but if we take the Green 
Bank measurement in Table 1 as an upper limit on the HVR V(H I), we obtain [0/H]hvr 
> —1.24 (the ionization correction should have little effect on this [0/H] estimate, as shown 
in Figures 11 and 13). This is consistent with the absolute metallicity derived for Complex 
C. However, the HV ridge contains substantially more highly ionized gas relative to the 
lower ionization stages. At the gas temperature and/or ionization parameter implied by 
the low-ion ratios in the HVR, Si IV and C IV have very small ion fractions and should 
be undetectable (see Table 5 in Sutherland & Dopita 1993 and Figure 12). Consequently, 
we conclude that the signiflcant high ion absorption lines observed in the high- velocity ridge 
originate in a separate phase from the low ionization stages. The similar component structure 
of the low and high ions in the HVR (see Figure 7 and Figure 8e-f) suggests that there is 
some association between the low-ionization and high-ionization phases; this could occur if 
the high ions arise in an interface between the low ion-bearing gas and a hotter ambient 
medium. 



-37- 



5.4. IVC Complex C/K 

The absorption-line measurements for the IVC C/K are substantially more uncertain 
than the HVC measurements due to line saturation and blending with adjacent components. 
Nevertheless, the 3C 351 spectrum indicates that C/K has a higher metallicity than Complex 
C proper and the HV ridge. The equivalent widths and column densities are significantly 
higher in the C/K component than in the high- velocity components (compare Tables 2 and 
3) despite the fact that the H I column is lower in C/K than in Complex C. Adopting 
the O I column from profile fitting (which provides better compensation for blending and 
saturation) and taking A^(H I) < 4 x 10^^ cm"^ (see § 2), we find that [0/H]c/k > -0.5. 
As we have shown in the previous sections, ionization corrections are likely to be important 
for the other species detected in C/K. We have modeled the ionization of the C/K gas using 
CLOUDY, and we find that the large deficit of N I with respect to Si II shown in Figure 9 
is predominantly due to ionization. At the ionization parameter that provides the best fit 
to the nominal O l/Si II ratio (i.e., log U — —3.6), the implied nitrogen underabundance is 
small and marginally significant, [N/0] = — 0.21q;3. The ionization-corrected abundances of 
the other detected species at this value of U are listed in the last column of Table 3. All 
abundance estimates in Table 3 are listed as lower limits since we only have an upper limit 
on Ar(Hl) (§ 2). 

6. Discussion 

The measurements and modeling presented in the previous sections have some inter- 
esting implications. We begin our discussion with some comments on the structure and 
confinement of Complex C (§6.1). We then compare the 3C 351 sight line to other nearby 
sight lines through this HVC (§6.2). We examine the case for metallicity variations in Com- 
plex C (§6.2.1), and then present evidence that the lower latitude region of Complex C is 
interacting more vigorously with the ambient medium (§6.2.2). Finally, we investigate the 
relationship between intermediate- and high- velocity gas observed toward H182 1+643 and 
the HVCs in the 3C 351 spectrum (§6.2.3). 

6.1. The Structure and Confinement of Complex C 

Our ionization models provide constraints on the physical conditions and dimensions of 
the absorber. The size of the absorber along the line-of-sight is L = N-n/nu- For the high 
(low) values of iV(0 I) in Complex C from Table 2, the CLOUDY models that provide the 
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best fits to the Complex C column densities (see § 5) have nu — 9 x 10~^ (3 x 10~^) cm~^, 
and L = 0.30 (2.0) kpc, and T f» 9300 (10, 700) K. These densities and temperatures suggest 
that the cloud may not be gravitationally confined. Using eqn. 3 from Schaye (2001), we 
see that a sclf-gravitating cloud with these Uh and T values would have a much larger size, 
L = 5.3 (10) kpc. The self-gravitating cloud size may change by a small amount depending 
on the absorber geometry, as noted by Schaye. However, even with this small scale factor, 
the uncertainties in the parameters derived from the CLOUDY calculations do not appear to 
be sufficient to reconcile the large difference between the expected size for a self-gravitating 
cloud and the size implied by the ionization models. 

The HVC could still be self-gravitating if it is predominantly composed of dark mat- 
ter. Following Schaye (2001), we have assumed that the fraction of the mass in gas fg ^ 
^b/^m — 0.16 to calculate the self-gravitating cloud size above. The cloud could be made 
self-gravitating by reducing fg to the order of 10~^ — 10~^. This may be unhkely, but it is 
not inconceivable. As an HVC plunges into a galaxy halo, ram pressure can separate the 
gas from its original dark matter halo (e.g., Quilis & Moore 2001). Wc may be viewing the 
small amount of residual gas in a dark matter halo that has already had most of its baryons 
stripped away. However, Quilis & Moore find that this requires the ambient medium to have 
a relatively high density, > 10~^ cm~^, which they deem "unrealistic for a Galactic halo 
component" . However, their analysis mainly considers HVCs at distances of > 100 kpc. If 
Complex C is only ~10 kpc away, the ambient halo density could be much higher, making 
this hypothesis more plausible. We note that the requirement of a small value for fg to make 
the cloud self-gravitating does not strictly require dark matter. However, searches for other 
components in HVCs such as stars (e.g., Willman ct al. 2002) or molecular hydrogen (e.g., 
Richter et al. 2002) have not produced detections despite good sensitivity. 

On the other hand, we may be observing the other side of this process, i.e., our gas 
may have been recently stripped out of a dark matter halo and is now largely devoid of dark 
matter. In this case, there is no particular need to confine the gas; the stripped gas could 
then be an ephemeral cloud that will rapidly evaporate in the hot halo. 

The purely coUisionally ionized model (§5.2.1) does not provide a direct constraint on the 
number density and hence the size of the absorber, but wc argue that the situation is similar. 
The coUisional model does provide an estimate of the H I ion fraction /h i. Substituting 

n„ ^ ^ (9) 

into Shaye's eqn. 3, we obtain 



T \ f 7V(HI) V' / /. 
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The best-fit from the coUisional ionization model with the high A^(0 I) has /h i = 0.43 and 
T = 10^-2 K, so with N{Yl I) = 4.2 x 10^^ we derive L = 130 kpc and = 2.5 x 10"^ cm'^ 
from equations 9 and 10 for a gravitationally confined cloud. The angular extent of Complex 
C (roughly 20° x 90°) implies that the transverse size of the entire complex is ~ 3.5 x 20 kpc 
if the cloud is 10 kpc away. Therefore L = 130 kpc requires an unlikely geometry: the cloud 
must be vastly larger along the line-of-sight than in the transverse direction. We reach the 
same conclusion with the parameters resulting from the lower O I column. This discrepancy 
could be reconciled, once again, by invoking a relatively small fraction of the mass in gas, fg. 
The implied geometry would also be more plausible if Complex C was considerably farther 
than 10 kpc, but there are arguments against placing this HVC much farther away (Wakker 
et al. 1999a; Blitz et al. 1999). We conclude that the coUisional ionization equilibrium 
model faces the same problems as the CLOUDY model. 

Of course, there are alternatives to gravitational confinement, and there is no require- 
ment for HVCs to contain dark matter at all. Magnetic confinement may be important 
(Konz, Briins, & Birk 2002). Pressure confinement by a hotter external medium is also a 
reasonable alternative. The widespread detection of high-velocity O VI absorption (Sem- 
bach et al. 2002) suggests that the Milky Way has a large, hot corona. Such a corona 
could provide the pressure confinement that we require. The densities and temperatures 
from the best CLOUDY models imply that the gas pressure in the Complex C component is 
p/k ^ 60 — 170 cm~^ K. If the external confining medium has Text ^ 10^ K, then its density 
^ext ^ 6 X 10~^ — 2 X 10~^ cm"^ in order to pressure confine the HVC. This density upper 
hmit is in agreement with density constraints from other arguments (e.g., Wang 1992; Moore 
& Davis 1994; Weiner & WiUiams 1996: Murah 2000; Briins, Kerp, & Pagels 2001). We note 
that these pressures are also consistent with constraints derived by Sternberg, McKee, & 
Wolfire (2002) on the pressure in the intragroup medium in the Local Group (based on H I 
structure properties in nearby dwarf galaxies). 

We can apply the same analysis to the high- velocity ridge component, but we only obtain 
limits since we only have an upper limit on A^(H I). For this high- velocity component, the 
best-fitting CLOUDY model implies an absorber thickness L < 2.6 kpc (decreasing A^(H I) 
also decreases L). With nn ~ 2.5x10"^ cm~'^ and T 11,300 K, the estimated size for a 
self-gravitating cloud (assuming fg — 0.16 again) is L = 11 kpc. In this case, decreasing fg 
by a factor of 2 — 3 reconciles the self-gravitating cloud size with the size implied by the 
ionization model. Gravitational confinement appears to be viable for the HVR, as long as 
the actual A'"(H I) is not too much lower than our upper hmit (see § 2). 

Summarizing this section, our constraints on the physical conditions and dimensions of 
the Complex C gas toward 3C 351 indicate that it is unlikely that this absorption arises in a 
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gravitationally confined cloud. It is much more probable that the gas is pressure confined by 
an external medium or is an ephemeral entity that will soon dissapate. The HVR could be 
gravitationally confined, but if the HVR and Complex C proper arc indeed closely related, 
then the physical processes affecting Complex C arc likely to affect the HVR as well, and 
the situation with the HVR may be similar to that of Complex C proper. 



6.2. Sight-Line Compcirisons 

6.2.1. Mrk 279, Mrk 817, and PG1259+593 

We can gain additional insights by comparing the 3C 351 results to other sight lines 
through the Complex C high-velocity cloud. Richter et al. (2001) and Collins et al. (2002) 
have published accurate column densities, including a variety of species, for several extra- 
galactic sight lines through Complex C. These papers have emphasized abundance results, 
but the observations have implications regarding the physical conditions/structure of the 
HVC as well. From these papers, the spectra of Mrk 279, Mrk 817, and PG1259-F593 have 
the highest S/N and provide the best constraints, so we shall concentrate on these sight 
lines. 

An important difference between 3C 351 and Mrk279/Mrk817/PG1259-F593 is that the 
latter sight lines have substantially higher H I column densities. The increased self-shielding 
resulting from the higher A^(H I) has advantages and disadvantages: on the one hand, 
the ionization corrections are smaller, but the ion ratios also provide less leverage on the 
density and physical conditions of the gas. Figure 14 shows an example. This figure shows a 
CLOUDY calculation, analogous to the model presented in §5.2.2, of metal columns vs. log 
U and log nn for the sight fine to PG 1259-1-593. The only difference from the 3C 351 model is 
that the H I column is higher (in accord with 21cm observations), and the overall metallicity 
has been reduced by a small amount to provide the best fit to the observed column densities. 
Because of the higher A^(H I) and greater self- shielding, most of the curves in Figure 14 are 
relatively flat; /(X*) = 1 over a large range of U for most species shown. Consequently, the 
derived abundances are insensitive to U, but also most of the measurements allow a large 
range of nn, L, and T. However, Ar I is a notable exception. As discussed by Sofia & Jenkins 
(1998), Ar I is susceptible to photoionization, and this can be useful for investigation of gas 
ionization. This can be seen in Figure 14. While most of the curves change by 0.1 dex or 
less, A^(Ar I) decreases by ~ 0.8 dex over the range of U shown in the figure. Unfortunately, 
the argon lines are relatively weak, and of the eight sight lines presented by Collins et al. 
(2002), Ar I is only detected toward PG1259+593. 
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Fig. 14. — Model of gas photoionized by the extragalactic UV background, as in Figure 12, 
but with parameters appropriate for the sight line to PG1259+593: log iV(H 1) = 19.92, 
Z = O.I2Z0, and solar relative abundances. The large points show the observed Complex 
C column densities in the direction of PG1259+593 from Collins et al. (2003). At log 
U = —3.9, all of the observed columns are within la of the model column densities except 
NX. 
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Taking advantage of the small ionization corrections, we confirm the main abundance 
results reported by Collins et al. (2002), with some additional important comments: 

1. The Complex C oxygen abundances^^ that we derive for these three directions show 
some dispersion: 

[0/H]Mrk279 = -0.591°:^^, (11) 
[0/H]Mrk817 = -0.50i°i^, (12) 

and 

[0/H]pGl259+593 = -OMt'o^ (13) 

However, including the 3C 351 oxygen abundance, we find that the weighted mean 
metaUicity for Complex C is <[0/H]> = -0.68 with a reduced xl = 0.85 (i/ = 3). 
Therefore the current measurements are consistent with a constant metallicity through- 
out Complex C with Z ^ 0.2Zq. Furthermore, the error bars in 11 - 13 do not fully 
reflect the A^(H I) uncertainties. For example, the H I 21 cm profile toward Mrk279 
is very complex, and it is difficult to separate the Complex CHI emission from the 
lower-velocity emission. Also, the MrkSl? A^(H I) was derived from the larger-beam 
LDS data, which introduces a systematic error of 0.3 — 0.5 dex. With these additional 
uncertainties, we conclude that there is currently no compelling evidence of oxygen 
abundance variations in Complex C. We cannot rule out the possibility that the oxy- 
gen abundances are spatially variable in Complex C, but better measurements are 
required to support this claim. 

2. Nitrogen is under abundant in Complex C: 

[N/H]Mrk279 < -1-1, (14) 

[N/H]Mrk817 < -1.2, (15) 

[N/H]Mrk876 = -l-20l°:l^, (16) 

and 

[N/H]pGl259+593 = -l-SSt'^oH (17) 

where the upper limits are at the Act level. Again, this indicates that intermediate- mass 
stars have not contributed significantly to the gas enrichment. The nitrogen detections 



-•^^Our oxygen abundances differ from those reported in Collins ct al. (2002) because we adopt the revised 
solar oxygen abundance reported by Holweger (2001) while Collins et al. adhere to the previous (O/H)© 
from Grevesse & Sauval (1998). We also include the uncertainties in the solar reference abundances reported 
by Holweger in the overall error bars. 
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toward Mrk876 and PG 1259+593 provide stronger evidence of abundance variations 
in Complex C than the oxygen measurements: these two sight hues yield a weighted 
mean of <[N/H]> = —1.40 with a reduced Xu ~ 5.60. Some caution is warranted, 
though, because N{N I) toward Mrk876 is entirely based on a weak and blended line 
(see Figure 8 in Collins et al. 2002), and the measurement may be more uncertain 
than the formal error bars suggest. 

3. The current error bars are too large to provide significant evidence of «— element 
overabundances. The model shown in Figure 14 provides a satisfactory match of the 
observed column densities with solar relative abundances, for example. The fact that 
the S II and Fc II column densities arc comparable suggests that the a— elements may 
be overabundant, but within current uncertainties, both [a/Fe] = and [ct/Fe] = +0.3 
provide acceptable fits to the observations. 

4. However, the A^(Fe II)/A^(S II) ratios do indicate that there is little depletion by dust 
in Complex C since dust strongly reduces the gas-abundance of Fe but has little effect 
on S in a variety of environments (e.g., Jenkins 1987; Sembach & Savage 1996). 

As noted above, most metals detected along these higher iV(H I) sight fines allow a wide 
range of densities and sizes. However, Ar I is useful. Toward PG1259+593, the best fit to 
the metals in Figure 14, including Ar I but excluding N I, has nu = 2.5 x 10"^, T = 8800 
K, and L — Nn/nu — 13.9 kpc. This impfied size is vastly larger than the size derived from 
the 3C 351 data. Furthermore, in the direction of PG1259+593, gravitational confinement 
appears to be quite viable: this only requires fg = 0.32, i.e., a factor of two larger than 
Schayc's fiducial value of 0.16. The upper limits on Complex C A^(Ar I) for the Mrk 279 and 
Mrk 817 sight lines do provide the following constraints. For Mrk 279, nn ^ 5 x 10~^ cm~^ 
and L > 2.5 kpc. For Mrk 817, nn < 1 x 10"^ cm"^ and L > 14 kpc. 



6.2.2. 3C 351: at the Leading Edge of Complex C 

Why do the sight lines to 3C 351 and PG1259+593 have such starkly different impli- 
cations regarding the size and confinement of the absorber? We propose a simple answer: 



I'^The Mrk876 N I column has also been estimated from the N I A1134.17 line observed with FUSE (Murphy 
et al. 2000). However, the feature at the expected velocity of N I A1134.17 in Complex C is strongly blended 
with low- velocity Fe II Al 133.67, and in fact the line is predominantly due to Fe II. Murphy et al. (2000) 
have subtracted this Fc II line based on other Fe II transitions in the FUSE bandpass, but the resulting 
iV(N I) estimate is considerably uncertain. 
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Fig. 15. — Map of the O VI column density in Complex C proper. The xy-plane shows the 
A^(H I = 2 X 10^^ cm"^ contour for Complex C from Hulsbosch & Wakker (1988), plotted vs. 
Galactic longitude and latitude. The z-axis shows N{0 VI), in units of 10^^ cm~^, measured 
towards nine QSOs/AGNs, from Table 8 in Sembach et al. (2002). The column densities are 
shown with filled circles, with a line marking the sight line location within Complex C in the 
contour plot in the xy-plane. A^(0 VI is also projected onto the xz-plane. The location of 
the 3C 351 sight line is marked with a +; O VI cannot be measured in the 3C 351 spectrum 
because of an extragalactic Lyman limit absorber at ^abs = 0.221. 



-45 - 




Fig. 16. — Map of the O Vl/H I ratio in Complex C proper, as in Figure 15. The z-axis 
shows the logarithm of A'"(0 VI)/A'"(H I) xlO^. The O Vl/H I ratios are also projected onto 
the xz-plane, and the location of the 3C 351 sight line is marked with a +. 
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Complex C is probably interacting with the thick disk/lower halo of the Milky Way much 
more vigorously in the direction of 3C 351 than in the direction of PG1259+593. This would 
occur naturally if the 3C 351 pencil beam is near the leading edge of the HVC, and this would 
explain several observations. First, the 21cm observations show a steep gradient in A^(H I) 
in the vicinity of 3C 351 (see Figures 1-3). This likely reflects the transition from the mostly 
neutral inner region to the fully ionized periphery of the HVC. Such a transition is expected 
at the leading edge of a cloud moving through an ambient halo (e.g., Quihs & Moore 2001). 
Second, ram pressure could separate the baryons from the underlying dark matter (Quilis 
& Moore 2001) and thereby alleviate the confinement problem discussed above. In this case 
there is no compelling requirement to confine the gas; this edge of the cloud would be rapidly 
evaporating as it plunges toward the disk. Third, observations of nine sight lines through 
Complex C reported by Sembach ct al. (2002) show that the O VI column densities at 
the velocities of Complex C are highest near 3C 351 (several of these nearby sight lines are 
marked in Figures 1-4). These O VI observations arc summarized in Figure 15, which shows 
iV(0 VI) measured toward the nine extragalactic sources marked on a map of Complex C. 
Note that these are O VI columns in Complex C proper only; recall also that O VI cannot 
be measured in the 3C 351 spectrum due to a Lyman limit absorber at Zabs — 

0.221 that 

severely attenuates the FUSE 3C 351 spectrum below ~1117 A. As shown in Figure 16, an 
even more pronounced trend is evident in the A^(0 VI)/A^(H I) ratio in Complex C. This ratio 
increases steadily with decreasing longitude and latitude, and the lowest-longitude Complex 
C sight line (Mrk 501; see Table 8 in Sembach et al. 2002) has an A^(0 Vl)/iV(H I) ratio 
which is ~43 times larger than the ratio measured toward the highest-longitude sight line 
(PC 1259-1-593). This trend also suggests that the lower- longitude section of Complex C is 
more strongly interacting with the ambient medium; as the gas is ionized to a greater degree, 
A'"(H I) will decrease while A^(0 VI) increases. The lower-longitude region is also at lower 
latitude, so this result is not surprising: as the HVC approaches the plane of the Milky Way, 
it is becoming more fully ionized and is probably ablating and dissipating, and the trailing 
(higher-longitude) portion of the cloud has not yet entered the higher density region of the 
ambient medium where the interactions are more vigorous, or at least has not suffered the 
effects of ablation for the same duration as the lower latitude/longitude regions. 

6.2.3. H1821+643 

An alternative to the interpretation presented above is that gas fiowing up out of the disk 
(e.g.. Galactic fountain/chimney gas) is colliding with Complex C at lower longitudes and 
latitudes, and the resultant shock-heating is ionizing and evaporating the cloud. Gibson ct 
al. (2001) and Collins et al. (2002) have suggested that such interactions might be required 
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to explain the metallicity variations observed in the HVC (but see point 1 in §6.2.1). In 
this regard, it is interesting to compare the 3C 351 to the H182H-643 sight hne, which is 
near 3C 351 but closer to the plane (see Figures 1-4). The H1821+643 line-of-sight passes 
just outside the lower-latitude 21cm boundary of Complex C (see also Figure 1 in Wakker 
& van Woerden 1997) but through the region of high-velocity emission associated with the 
Outer Arm. A good spectrum of H182 1-1-643 could reveal gas flowing from the plane toward 
Complex C. 

Savage, Sembach, & Lu (1995) have discussed Galactic absorption features in the spec- 
trum of H1821-I-643 based on an intermediate-rcsohition GHRS spectrum. Subsequently, 
Tripp, Savage, & Jenkins (2000) obtained a higher resolution STIS spectrum of H1821-I-643 
with the same mode (E140M) used to observe 3C 351, so a detailed comparison of these 
sight lines can be made without confusion due to differing spectral resolution. Figure 17 
compares selected absorption profiles observed toward H1821-I-643 (upper profile in each 
panel) to the same line from the spectrum of 3C 351. For reference, the velocities of the 
HVR, C, and C/K components are marked at the top of each panel. Some of the profiles 
appear to be remarkably similar, e.g., the Fe II lines in Figure 17d. We defer a full analysis 
of the H1821-I-643 data to a later paper, but a few comments on this figure are worthwhile. 

The V = —80 km component: Complex C/K. Both sight lines show components 
near the velocity of IVC C/K with similar equivalent widths. Toward H1821-I-643, Savage 
et al. (1995) associate this absorption with the Perseus spiral arm, which has a similar 
velocity at b ^ and appears to extend ~ 1 kpc above the plane based on 21cm and 
Ha emission (Kepner 1970; Reynolds 1986). We have derived [0/H] > —0.5 for the C/K 
component toward 3C 351, implying a higher metallicity than the adjacent Complex C and 
HVR components. The higher C/K metallicity suggests that this gas may indeed have 
originated in the disk and is being driven into the halo. If the Perseus arm is ~ 3 kpc 
away and the 3C 351 C/K component is associated with Perseus, then this gas is ~ 2 kpc 
above the plane. This is consistent with the ^—heights that gas is expected to attain in 
some Galactic fountain models (e.g., Houck & Bregman 1990). Very recently, Otte, Dixon, 
& Sankrit (2003) have detected O VI emission at ^lsr = ~50 ± 30 km s~^ from / = 95.4°, 
b = 36.1°, a direction very close to the 3C 351 and H1821-I-643 sight lines. They also show 
bright Ha filaments extending up from the plane in this direction. They attribute both the 
O VI and Ha emission to an outfiow from the Perseus arm. Regardless of the nature of 
the gas, it is quite likely that the absorption lines a,t v ^ —80 km toward 3C 351 and 
H1821-I-643 are related to the emission reported by Otte et al. (2003). 

The V — —130 km component: Complex C vs. the Outer Arm. Strong absorption 
fines are present at the velocity of Complex C in the spectra of both H1821-I-643 and 3C 
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Fig. 17. — Comparison of normalized absorption profiles observed in the directions of 3C 
351 (lower profile in each panel) and H1821+643 (upper profile in each panel). The plots 
compare the (a) O I A1302.17, (b) Si II A1304.37, (c) Al II A1670.79, (d) Fe II A1608.45, (e) 
Si III A1206.50, and (f) C IV A1548.20 transitions. The H1821+643 spectra have been offset 
by +1 for clarity. These sight lines probe the low-latitude edge of Complex C (see Figures 1 
- 3 as well as Figure 1 in Wakker & van Wocrdcn 1997). For reference, the velocities of the 
high- velocity ridge (HVR), Complex C, and JVC C/K are marked at the top of each panel. 
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351. However, unlike IVC C/K, the line equivalent widths in the H1821+643 spectrum are 

substantially larger than the same lines toward 3C 351 at v = —130 km s~^ (see Figure 17). 
The H1821+643 and 3C 351 21 cm H I column densities at this velocity are similar: Wakker 
et al. (2001) report A^(H I) = 3.3 x 10^^ toward H1821+643, compared to iV(H I) = 4.2 x 10^^ 
toward 3C 351 (both measurements are based on Effelsberg observations). However, there 
is a gap in the 21 cm emission between the two sight lines (see Figures 1 and 3), and the 
relationship (if any) between the 3C 351 and H1821+643 HVCs at this velocity is unclear. We 
list in Table 4 the equivalent widths and column densities of absorption lines aiv — —130 km 
in the H1821+643 spectrum, measured as described in § 4. Several pixels in the core of 
the O I hue toward H1821+643 approach zero flux. This hne may be significantly saturated, 
so we consider the O I measurements highly uncertain. The most conservative constraint on 
A^(0 I) is a lower limit from direct integration, but we also provide an estimate from profile 
fitting in Table 4, which can correct somewhat for saturation but with large uncertainties. 
We also summarize the implied abundances in Table 4 as before, without and with ionization 
corrections applied (using the ionization model presented below). 

The larger equivalent widths toward H1821+643 at wlsr ~ —130 km s^^ ostensibly indi- 
cate that the metallicity at this velocity is greater in the H182H-643 component than in the 
3C 351 HVC since the H I columns are similar toward both QSOs. If we neglect ionization 
corrections, we indeed obtain high abundances from the column densities in Table 4, e.g., 
[O I/H I] > -0.6, [Si II/H I] = 0.13, [Al II/H I] = 0.0, and [Fe II/H I] = -0.08, but with a 
surprisingly low nitrogen abundance, [N l/H I] = —0.97. This would be an unusual abun- 
dance pattern; N underabundances are expected in low- metallicity gas, but as Z — > Zq, 
the contribution from intermediate-mass stars is expected to bring the N relative abundance 
more in line with the solar value. However, we have shown in § 5 that ionization corrections 
must not be neglected when A^(H I) is at the level observed toward H1821+643. When ion- 
ization corrections are applied, the observed columns can be reconciled with a lower overall 
metallicity and smaller N underabundances. Figure 18 shows a CLOUDY model in satisfac- 
tory agreement with the observed columns with Z — 0.24Zq and only small underabundance 
of N and Al. The model shown in Figure 18 works if A'"(0 I) is close to the value from direct 
integration, i.e., A'"(0 I) ^ 14.7. Higher O I column densities would require higher metal- 
licity, lower ionization parameters (to match, e.g., the O I/Si II ratio), and greater nitrogen 
underabundances, as can be seen from Figure 18. A slight depletion of Fe might also be 
required in higher metallicity models. Due to the large range of U allowed by the current 
data, we cannot usefully bracket the size of the absorber at vlsr = —130 km s~^. 

At this juncture, it remains possible that the Outer Arm high- velocity gas and Complex 
C have similar metallicities, but the data also allow the Outer Arm to have a substantially 
higher metallicity than Complex C if N is under abundant. Consequently, it is difficult to 
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Table 4. High- Velocity Absorption Lines toward H1821-I-643: The Outer Arm*^ 



Species Aq <v> Wx log A^a log iVpf [XVH I]^ [X/H]^ 

(A) (km s-i) (mA) 



I 


1302.17 


-131 ± 1 


183.6±5.4 


> 14.70^^ 


10.iO_Q 24 


> -0.55 


> -0.57 


N I 


1199.55 


-136 ±5 


40.0±9.3 


13.47±0.11 


13.40 ±0.09 


-0.97 ±0.17 


-0.73 ±0.17 


Si II 


1304.37 


-128 ± 1 


134.0±4.4 


14.18 ±0.02 




U.iC5_o.i2 






1526.71 


-129 ± 1 


196.6±6.0 


14.15 ±0.04 




^•^^•^-0.12 




Al II 


1670.79 


-128 ± 1 


222.8±10.0 


12.98±0.04 


13.02 ±0.04 


0.00 ±0.11 


-0.80 ±0.11 


Fe II 


1608.45 


-131 ±1 


79.3±5.4 


13.88±0.04 


13.89±0.04 


-0.08 ±0.12 


-0.60 ±0.12 



^See Table 2 for definitions of the quantities in this table. The equivalent widths and apparent column 
densities are integrated over the velocity range of the high-velocity component afiiliated with the Outer Arm, 

from wlsr = —169 to —99 km s~^. 

''Implied logarithmic abundance if ionization corrections are neglected, i.e., [X^/H I] = log A^(X*)/iV(H I) - 
log (X/H)0. 

'^Logarithmic abundance obtained by applying the ionization correction from the CLOUDY model shown 
in Figure 18 and discussed in § 6.2.3. Error bars include column density uncertainties and solar reference 
abundance uncertainties but do not reflect uncertainties in the ionization correction. In this case, the ionization 
correction allows a large range of abundances (see § 6.2.3). 

''Saturated absorption line. 

''Due to considerable saturation, the formal error bars from Voigt-profile fitting may underestimate the true 
uncertainty. 

^Simultaneous fit to the Si II A1304.37 and A1526.71 lines. 
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Table 5. High- Velocity Absorption Lines toward H1821+643: The -212 km s"^ 

Component^ 



Species 


Ao 

(A) 


<v> 
(km s~^) 


(ml) 


log /^a 


log Npf 


I 


1302.17 




3.8 ±4.7 


< 13.38'' 




C II 


1334.53 




7.6 ± 3.4 


< 12.84'' 




Si III ... . 


1206.50 




-3.9 ±8.6 


< 12.21^^ 




Si IV 


1393.76 




10.0 ±4.6 


< 12.32'' 




C IV 


1548.20 


-212 ± 1 


52.4 ±5.5 


13.22 ±0.05 


13.24 ±0.03^ 


N V 


1238.82 




7.4 ±5.4 


< 13.01'' 




VI 


1031.93 


-192<^ 




13.72±0.14'^ 





'^See Table 2 for definitions of the quantities in this table. The equivalent 
widths and apparent column densities are integrated from wlsr = —234 to 
-184 km s-^ 

''4ct upper limit derived from the equivalent width limit assuming the linear 
curve-of-growth applies. 

'^Column density from a joint fit to the C IV 1548.20 and 1550.78 A transi- 
tions. 

'^O VI line centroid and column density from Sembach et al. (2002). 
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Fig. 18. — Model of gas photoionized by the extragalactic UV background, as in Fig- 
ure 12, but with parameters appropriate for the Outer Arm component in the spectrum 
of H1821+643. The large points show the observed Outer Arm column densities toward 
H1821+643. At log U = —3.6, the model is in reasonable agreement with the observed 
column densities. However, A^(0 1) may be underestimated due to saturation; to fit the 
observations with a higher O 1 column density, U must decrease and the metallicity must 



mcrease. 
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draw firm conclusions about the relationship between the Outer Arm and Complex C. Weaker 
O I lines in the FUSE bandpass may help to more tightly constrain the metallicity of the 
Outer Arm and thereby clarify this relationship. 

The V — —200 km s"^ component: the high-velocity ridge. Finally, it is interesting 
to compare the 3C 351 and H1821+643 sight fines at the velocity of the HVR component, 
vlsr ~ —200 km s~^. While the HVR absorption is clearly detected in a variety of transitions 
toward 3C 351, most of the H1821+643 profiles show nothing at this velocity (the Si III 
profiles in panel (e) of Figure 17 show the contrast between the sight lines at this v). However, 
Savage et al. (1995) tentatively identified C IV absorption near this velocity in a GHRS 
spectrum of H1821+643. As shown in Figure 19, we confirm this detection: both lines of 
the C IV AA1548.2, 1550.8 doublet are well-detected at = -212 km s^^ These are the 
only lines that we detect in the STIS spectrum at this velocity. The C IV equivalent width 
and column density measurements are listed in Table 5 along with upper limits on selected 
dominant ions and higher stages. However, Oegerle et al. (2000) and Scmbach et al. (2000) 
detected O VI at this velocity in FUSE spectra of H1821+643. The O VI 1031.9 A fine is 
blended with low-velocity H2, but Sembach et al. (2002) have removed the H2 line (based 
on fits to other H2 fines) and find log N{0 VI) = 13.72 ± 0.14. 

These measurements have useful implications. The fact that iV(0 VI) ^ iV(C IV) 
toward H1821+643 would indicate that T > 10^ K if the gas were collisionally ionized and 
in equilibrium (see, e.g.. Figure 7 in Tripp & Savage 2000). However, this temperature is not 
compatible with the widths of the C IV lines. Figure 19 shows our Voigt-profilc fit to the HVR 
features in the H1821-I-643 spectrum; a single component with h = 9 ±2 provides an excellent 
fit. This implies that T < 10^-^^ K. The C IV line width is only marginally consistent with 
the O Vl/C IV ratio given the uncertainties in A^(0 VI) and b{C IV). It seems more likely 
that either (1) much of the O VI absorption is not associated with the C IV-bearing gas (in 
which case the O Vl-gas must be relatively hot to satisfy the lower limit on O VI/C IV), or 
(2) the gas is not in ionization equilibrium. The different centroids of the C IV and O VI lines 
(see Table 5) support the idea that at least some of the O VI originates in different gas. It is 
unlikely that the C IV + O VI gas is photoionized because of the long pathlengths required 
(e.g., Sembach et al. 2002), but it is possible that the gas is collisionally ionized but out 
of equilibrium because it is coofing faster than it can recombine (Edgar & Chevafier 1986). 
Heckman et al. (2002) have recently shown that a non-equilibrium, radiatively cooling gas 
model provides a good fit to the O VI column densities and 6— values in HYCs (as well as 
other contexts). Their model also predicts A^(0 VI) » A'"(C IV), as observed in the —200 
km s-^ feature toward H1821+643. 

Savage et al. (1995) have noted that i'lsr ~ —200 km s"^ is near the expected velocity 
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Fig. 19. — Continuum-normalized C IV absorption profiles observed toward H1821-I-643. 
The upper histogram shows the A1550.8 transition and the lower histogram is the A1548.2 
line. The fit to the component at i'lsr = —212 km s~^ is overplotted with a sohd line. 
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for distant Milky Way gas in a corotating disk/halo in the direction of H1821+643; the 
velocity would then imply a large Galactocentric distance (see their Figure 2). However, 
it is quite possible that this high-velocity feature is associated with the high-velocity ridge 
observed toward 3C 351 and other sight lines. In this case, the fact that the feature is only 
detected in C IV and O VI would indicate that the HVR has an extended ionized periphery. 
Evidence for highly-ionized layers on the surface of HVCs has been reported for other clouds 
(e.g., Sembach et al. 1999). We have argued that the absorption hnes toward 3C 351 indicate 
that the HVR and Complex C proper are related. If the ^lsr = —212 km s~^ component 
toward H1821-I-643 is also part of the HVR, then this gas is hkely much closer than imphed 
by its velocity and the assumption that it is corotating (various arguments suggest that 
Complex C is ~10 kpc from the Sun). 



7. Summary 

We have investigated the physical structure, conditions, metallicity, and nature of the 
high- velocity cloud Complex C using high-resolution recordings of absorption lines in several 
directions through the cloud. Our study is mainly based on STIS echelle spectroscopy of 3C 
351, which shows a wide variety of absorption lines from Complex C proper (vlsr = —128 
km s~^) as well as an intermediate- velocity cloud (C/K, t^LSR = ^82 km s~^) and a higher- 
velocity component that we refer to as the high- velocity ridge (^lsr = —190 km s~^). We 
also make use of other sight lines through Complex C from the literature as well as new STIS 
echelle observations of H1821-I-643, which is just outside the 21 cm boundary of Complex C 
but shows absorption lines at very similar velocities to those seen toward 3C 351. Prom our 
analysis, we reach the following conclusions: 

1. The high- velocity ridge is closely related to Complex C proper. This is suggested by 
the similar morphologies of the HVR and Complex C; the HVR has a similar shape and 
is roughly centered on the larger, lower- velocity Complex C. This idea is supported by 
the absorption lines, which show remarkably similar column density ratios in the HVR 
and Complex C. However, very little high-ion absorption is detected at the velocity of 
Complex C proper, but Si III, Si IV, and C IV absorption is strong and easily detected 
in the HVR. Moreover, the high-ion absorption lines in the HVR have centroids and 
hne widths that are similar to those of the low ions in the HVR. We show that the high 
and low ionization stages in the HVR cannot arise in the same phase, and we conclude 
that the high ion absorption occurs in an interface between the low-ionization phase 
and a hotter ambient medium. 

2. The relative abundances in all components with A'"(H I) < 10^^ cm~^ indicate that 
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ionization corrections are important. We consider coUisional ionization equilibrium as 
well as CLOUDY photoionization models to derive ionization corrections, and we find 
in both cases that Z = (0.2 ±0.1)^0 in Complex C in the direction of 3C 351. We also 
find that nitrogen must be underabundant toward 3C 351. The low metallicity and 
nitrogen underabundance indicate that Complex C is not ejecta generated in a Galactic 
fountain. It seems more likely that Complex C is either tidally stripped material from 
a satellite galaxy (analogous to the Magellanic Stream), or that it is gas with a more 
distant extragalactic origin. The absolute metaUicity of the 3C 351 HVR component 
is less constrained because we only have an upper limit on the HVR H I column. The 
lower limits on the HVR metallicity are consistent with the metallicity derived for 
Complex C proper. 

3. We find similar oxygen abundances at Complex C velocities toward Mrk 279, Mrk 
817, and PG1259+593. While there is some dispersion in the [0/H] measurements, 
within the current uncertainties the measurements are fully consistent with a constant 
metallicity throughout the HVC. These sight lines as well as the Mrk 876 sight line 
also provide strong evidence of N under abundances in Complex C. Comparison of 
the nitrogen abundances toward Mrk 876 and PG1259+593 provides the strongest 
evidence of spatial abundance variability in Complex C, but this requires confirmation 
with additional measurements. There are hints of a— element overabundances in some 
directions, but the evidence is not statistically significant. 

4. The derived iron abundances indicate that Complex C contains little or no dust. 

5. Toward 3C 351, the Complex C absorber size, density, and temperature implied by the 
ionization models indicate that the gas is not gravitationally confined, while toward 
other sight fines such as PG1259+593, the implied path length through Complex C 
is much larger and gravitational confinement is viable. Pressure confinement by an 
external medium may play a role. However, we suggest that the 3C 351 sight line passes 
through the leading edge of Complex C, which has been ablating and dissipating as it 
approaches the plane of the Milky Way. This idea is supported by O VI observations, 
which show the highest O VI column densities near the 3C 351 sight fine. Furthermore, 
the O Vl/H I ratio increases dramatically with decreasing longitude and latitude within 
Complex C, which also suggests that the lower longitude and latitude regions are 
interacting more vigorously with the ambient medium. 

6. To look for evidence of outfiowing Milky Way gas that might interact with Complex C, 
we compare the STIS echelle observations of the sight fines to 3C 351 and H1821+643, 
which are both near the lower-latitude edge of Complex C at Z 90°. We find that 
the intermediate-velocity gas observed toward both of these QSOs likely has a higher 
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metallicity and may indeed be a fountain/chimney outflow from the Perseus spiral arm. 

Unfortunately, the results are less clear for the high-velocity gas toward H1821+643. 
The H1821-I-643 HVC could have a similar metallicity to that derived for Complex C, 
or it could have a much higher abundance depending on the O I column and ionization 
correction, both of which are highly uncertain. 

As usual, we are grateful to Gary Ferland and collaborators for the development and 
maintenance of CLOUDY. The observations reported here were obtained by the STIS Inves- 
tigation Definition Team, and this research was supported by NASA through funding for the 
STIS Team, including NASA contract NAS5-30110. TMT appreciates additional support for 
this work from NASA Long Term Space Astrophysics grant NAG5-11136 as well as NASA 
grant GO-08695.01-A from the Space Telescope Science Institute. BPW was supported by 
NASA grant NAG5-9024. 
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